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Abstract 

Using a connection between the q-oscillator algebra and the coefficients of the 
high temperature expansion of the frustrated Gaussian spin model, we derive an 
exact formula for the number of closed random walks of given length and area, on 
a hypercubic lattice, in the limit of infinite number of dimensions. The formula is 
investigated in detail, and asymptotic behaviours are evaluated. The area distri- 
bution in the limit of long loops is computed. As a byproduct, we obtain also an 
infinite set of new, nontrivial identities. 
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1 Introduction 



Let us consider on an hypercubic lattice of dimension D all closed paths of length 2k, 
starting from a given vertex. An interesting and still unanswered question concerns the 
counting of such paths according to the area they enclose. Besides being a challenging 
combinatorial problem on its own right, the question has some relevance also from a 
physical point of view: considering for instance the hopping of a charged particle around 
a two dimensional lattice in a uniform magnetic field |l], ||, || , its energy spectrum displays 
an intricate hierarchic structure || , whose properties are closely related to the area versus 
length distribution of closed paths p. Such distribution finds applications also in polymer 
physics H 0, and in explaining anomalous magnetoconductance ||. A different but 
somehow related problem is the one of counting random walks according to the number 
of distinct sites they visit P, |10[|: relatively compact loops are expected to enclose a small 
area and to visit relatively few sites, and viceversa. This counting problem is relevant 
when investigating random walks on lattices with static traps. 

Coming back to the initial question of counting loops according to their enclosed 
area, in the case of a two-dimensional lattice the answer was given only in the limit of 
asymptotically large k in [[11]], and the first subleading correction was provided in | I2[ . In 
higher dimension, the problem is in general even more difficult. Here we want to address 
the question in the limit of infinite dimensionality of the lattice, a very peculiar situation, 
in which however some simplifications occur, and an exact answer can be given. The latter 
is to be interpreted as a "mean-field" approximation to the exact counting on a lattice 
of finite dimensionality D. Such a formulation of the problem was first considered by 
Parisi et al. |TB|, [T3j]. They were investigating spin models with frustration but without 
any quenched disorder, in order to test the conjecture that such deterministic models 
could behave at low temperature as some suitably chosen spin-glass model with quenched 
disorder. They considered the frustrated Spherical and XY spin models in the limit of 
large dimensionality D of the lattice, where the saddle point approximation becomes 
exact. In their analysis of these models, they showed how the high temperature expansion 
(i.e. loop expansion) can be nicely rewritten by using the g-oscillator algebra \TE\, where 
q measures the frustration per plaquette. Similar g-deformed algebraic relations have also 
appeared in the Hofstadter problem of quantum particles hopping on a two-dimensional 
lattice in a uniform magnetic field JTBJ. This problem is closely related to the frustrated 
spin models, which can be considered as simplified models of hopping in the large D 
and classical limits. The use of somehow related non-commutative geometry techniques 
appeared ClS db crucial ingredient also in Ref. |J12|. 



The problem we address and solve in this paper may be stated as follows: given an 
infinite dimensional hypercubic lattice, what is the number Gf. i of loops of length 2k, 
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starting from a given vertex, enclosing a minimal area of exactly I plaquettesQ? 

The paper is organized as follows. In next Section we summarize some results of 



14j] which are then used in Section 3 to obtain an exact expression for the counting 



numbers Gk,i- In Section 4 a generating function for these counting numbers is presented. 
In Section 5 a probability distribution is naturally associated to the Gk/s, and a sys- 
tematic procedure is built to compute its moments of arbitrary order. The asymptotic 
behaviour of this probability distribution for large k is evaluated in Section 6, while Sec- 
tion 7 is devoted to the presentation of some new, nontrivial identities. The last Section 
contains our conclusions, while some technical details are relegated into the Appendix. 



2 Frustrated Lattice and ^-oscillator Algebra 

In Refs. |P~3| , 13], Parisi et al. considered the frustrated Gaussian, Spherical and XY 
models on a D-dimensional hypercubic lattice, in the large D limit. Among other things, 
they unveiled a remarkable connection between the coefficients of the high temperature 
expansion (i.e., loop expansion) and the g-oscillator algebra |L5|, where q measures the 



frustration per plaquette and varies continuously on the real interval [—1,-1-1], between 
the fully-frustrated case (q = —1, fermionic algebra) and the ferromagnetic case (q = 1, 
bosonic algebra). For our purposes, the relevant part of the Hamiltonians describing the 
models considered in Refs. |13|, [L4j is: 

H = —m E^*** + h - c - • C 1 ) 

The complex field <ftj £ C is defined on the sites (labelled by j) of a D- dimensional hyper- 
cubic lattice; frustration is induced on the lattice through the nearest-neighbour couplings 
Ujk, which are complex numbers of modulus one and satisfy the relation Ujk = U^; they 
are the link variables of a background Abelian lattice gauge field, chosen in such a way 
as to produce a static and constant magnetic field, suitably oriented to give the same 
magnetic flux ±B for any plaquette of the lattice (the product of the four U's around 
the plaquette being e ±lB ). Such a magnetic field having the same projection over all 
the axes, up to the sign, in Refs. [|T^, 13 these signs were chosen randomly in order to 



avoid the selection of a preferred direction. The usual unfrustrated ferromagnetic spin 
interaction is obtained for B = 0. Non-vanishing values of B induce a frustration around 
each plaquette, which is maximal for B = n, the fully-frustrated case. 



^When facing this problem of area versus length distribution of closed random walks, the area is most 
often considered in its algebraic sense, i.e. as oriented; in our case, however, an average over orientations 
is intrinsecally performed in the approach (see Section 2), and we are therefore referring to the absolute 
value of the area. 
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In the framework of the high temperature expansion, the free energy of such models 
is expressed as a sum over the contribution of loops of increasing length 2k: 

oo o2k 

P F = J2 2k Gk ■ ( 2 ) 



Each loop encloses a number of plaquettes; in the case of the models considered in [13, 14 



for each loop the magnetic field yields a weight, proportional to exp(iBA), where A is 
the sum of plaquettes with signs depending on the orientations. The total contribution 
of all loops of length 2k is given by Gk- Due to the average over orientations and loops, 
the quantity Gk is a polynomial in the variable q = cosB, the coefficient Gk,i of q l being 
given by the number of loops of length 2k and area / 0. The infinite dimensionality of 
the lattice ensures that for any given loop, no two plaquettes of the subtended minimal 
surface will lie on the same plane; this in turn accounts for such a simple averaging over 
orientations, performed indipendently for each plane. The order of polynomial Gk{q) is 
given by kik — l)/2, the maximal area encloseable by a loop of length 2k. The coefficients 
Gk,i can also be interpreted as the number of Feynman diagrams with 2k external points, 
which are joined pairwise by lines (propagators) intersecting I times. These diagrams also 
occur in the topological (large N) expansion of Matrix Models ||7j , where the planar limit 
corresponds to no intersections, i.e. to the q = case. Equivalently, in simple graphical 
terms, they just can be seen as the number of way of connecting pairwise 2k points on 
a circle with k segments intersecting exactly / times. In the following we shall refer to 
this last picture^ In Ref. [I3|, the enumeration of such diagrams was investigated. In 
particular, a recursion relation was found for the coefficients of the polynomial Gk(q) - 
a sort of Wick theorem - which can be nicely expressed in terms of the algebra of the 
g-oscillators a q , a q : 

a q a\ - qa\a q = 1 . (3) 

These operators [15j act on the Hilbert space spanned by the vectors: |m), m = 0, 1, 
as follows, 



a q \m) = y/[m] q \m - 1) , a q \0) = 



41 



1 -q r 



%\ m ) = y[m + l] q \m+l) , [m] g = - _ - . (4) 

Using the recursion relation, the weighted multiplicities of the diagrams of eq. (fj) were 
neatly written as an expectation value over the ground state of the g-oscillators |I3L |i~4"f: 



G k (q) = (0\(al + a q ) 2k \0) . (5) 



■■■The area of a loop is defined as the minimal area of a surface of lattice plaquettes having that loop 
as boundary. 

§This formulation of the problem was considered e.g. by Linus Pauling ]l8| when trying to simplify 
the calculation of matrix elements involved in Slater's treatment of the electronic structure of molecules. 
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The authors of Refs. |L3], [L4| did not exploit the consequences of this result, preferring 
to turn themselves to numerical investigation of the models they were interested in. As a 
consequence, till now only the two limiting cases Gfc(O) and Gfc(l) were explicitly known; 
when q — 1 it is only a matter of counting the way of connecting 2k points on a circle, 
with no restrictions, and this is simply the number of pairings of 2k objects: (2k — 1)!! . 
When q = we are in fact evaluating the planar limit of the zero-dimensional 2/c-point 



Green function of a Matrix Model in the limit of vanishing interaction |17j. In simple 
graphical terms this correspond to the number of way of joining pairwise 2k points on a 
circle with no intersection. In other words, this is just one of the many possible definition 
of Catalan numbers (see, e.g. |[9|]), given as: 

But the number Gkj of way of connecting pairwise 2k points on a circle, with exactly 
/ intersections (i.e. the coefficient of q l in Gk(q)) is not so easily accessible. In Refs. 
[ PI jjj , they were found by direct enumeration of the graphs on a computer. In pOf . 



a generating function for coefficients Gk,i was obtained,proposed, but then the explicit 
evaluation of such coefficients relied upon heavy symbolic manipulations. 

3 Exact Solution of the Enumeration Problem 

We shall now present a simple, easy to evaluate, formula for a generic coefficient Gk,i- 
To this purpose let us first introduce the x q coordinate representation, x q = a) q + a q , 
x q \x) = x\x), which is given by the so-called continuous g-Hermite polynomials [^2. 23 . 
These are defined by: 

H n (x) = (x\n) C n , C n = ([n] ? !) 1/2 C , (7) 
where the normalization constant Cq is fixed by H (x) = 1 and the g-factorial is 

M, ! = H [n-l) q ... [1], , [1], = [0], = 1 . (8) 
These polynomials satisfy, of course, a three-term recursion relation in the index n: 

xH n (x) = H n+l (x) + [n) q H n _ x (x) n > 1 . (9) 

and a convenient convenient 



x ranges over the interval x G —2/^/1 — q, 2/\/l — q 
parametrisation is 



x= 2 cosfl, 0G [0,71-1. (10) 

V 1 - Q 
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More properties of these g-Hermite polynomials can be found in Ref. |[22|| , where they 
are defined as H n (cos9) = (1 — q) n ^ 2 H n (x). The most important property for us is the 
orthogonalizing measure v q {pc) |^Tj, 0, [2! 



v q {x) dx H n (x) H m (x) = <5 n , m [n) q \ , (11) 



Z7T \ 71 



T 



n n=0 



where Qi(z,t) is the first Jacobi theta function 24 



0! ( z , r) = e i7rr(n+ 5) 2+27 ™(™-4)( 2 -§) , q = e 27riT . (13) 

neZ 

It is worth emphasizing that, interestingly enough, measure ( |I~2| ) interpolates continuously 
between Wigner semi-circle law for Gaussian Matrix Models, and Gaussian distribution, 
as q varies from to 1~; this reminds us of the behaviour of the density distribution of 
eigenvalues for Gaussian Ensemble of Random Matrices, as the dimension of the matrices, 
N, varies from oo to 1. 

Polynomial Gk{q) may now be rewritten as 

G k {q) = C' ^ v q (x) x 2k dx . (14) 



Indeed, using explicit form ( |12]) for the integration measure, performing the integral and 
playing a little bit with indices, we may write 

cm - (^) g(-D' |ii ^ • m 

We now only need to perform explicitly the division of the polynomial of degree 
kik + l)/2 defined by the sum in the previous formula. To this purpose, let us state the 
following 

Theorem^]: Let P(q) = Y2i=oPiQ l be an integer coefficient polynomial of degree N 
in q, exactly divisible by (1 — q) k . Then A(q) = P(q)/(1 — q) k is an integer coefficient 
polynomial of degree N — k whose coefficients are simply expressed in terms of the pi 's as 
follows: 

A(q)=Ya iq l , ai = E(* + i~! -< V (16) 

i=o i=o V K 1 / 



^We are unfortunately unable to quote any reference. We are of course convinced this theorem has 
been enounced long time ago, indeed it is just a corollary of Ruffini's rule. The only proof we have been 
able to give, absolutely inelegant and unworthy of being presented here, is by direct verification. 
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We therefore readily get a simple closed expression for Gk,f- 

2i+l 



G 



k,i 



where L 



EC- 1 ) 4 



i=0 



fc + Z- 1 - i(i + l)/2 
fc — 1 




/ < 



k(k-l) 



(17) 

1 + yl + 8/)/2 is the largest integer i satisfying %{% + l)/2 < /. 

J i.p. 

Because of their definition as counting numbers, all integers G% i should be positive. This 
fact is indeed not apparent from fll7l) , and we are not able to present a rigorous proof of 
this statement, nevertheless we are convinced of its validity. Sensibleness arguments rely 
upon the intrinsic positivity of (|5|) and its derivatives with respect to q, for any value of 
q in the interval [0, 1). 

We have computed explicitly with Mathematica P5], through formula fllTD , which is 
very efficient, all values Gk,i for k ranging from 1 to 9, and verified that they indeed 



match the results of Ref. [13 which were found by direct enumeration of the graphs on 
a computer. 



4 Generating Function 

Let us define a generating function for the combinatorial numbers Gkf- 

oo oo 

R{z, q) = J2 G ^) zk = E G M ■ (18) 

k=0 k,l=0 

The running of index I may or may not be restricted to the range < I < k(k — l)/2, 
without loss of generality: coefficients Gk,i automatically vanish for I > k(k — l)/2. From 
expression (|T^) , we may readily write a generating functions for coefficients Gk{q)'- 

/•2/vT=« 1 

R{z,q) = / v q (x) =■ dx . (19) 

Inserting for the measure its explicit expression (0), and performing the integral, we 
obtain: 

z 



I -i _ oo / \ 2n+l 

R(z,q) = J!—± E ("I)" <f (n+1)/2 



» n=0 

where we introduced 



7\ 



(20) 



1 - VI -4a 

^ («) = 7TT • 21 

For q = function (^Dj) reduces to the standard generating function of Catalan numbers. 
The singularities of R(z, q) in the complex z-plane are completely determined by those of 
the simpler function K(zf{\ — q)). Actually the series is very well convergent for |g| < 1 
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as any Jacobi theta function. Moreover, for \q\ = 1, it is a geometric series, which is still 
convergent because |i^(a;)| < 1 for \z\ < z c = (1 — g)/4. 

Generating function (p0|), as shown in Ref. |2(| in a different context, can be inter- 
preted physically as the internal energy at temperature T = l/\/z of a Gaussian spin 
model defined on the sites of the infinite-dimensional frustrated lattice introduced in Ref. 

0,111. 



5 Moments of the Distribution 



The interpretation of coefficients G k j naturally suggests to introduce, for each positive 
integer value of index k, the following normalized (discrete) probability distribution: 



Pk{l) 



Gk 



G 



k,i 



G 



k.i 



fc(fc-l)/2 



=0 



G 



k.l 



G k {\) (2k-l)\\ 



lb = 0,1,... 

Z = 0, 1, . . . , — l)/2 



(22) 

whose interpretation is obvious: given an arbitrary diagram with 2k legs, the probability 
for it to have exactly / crossings is simply given by P k (l); equivalently, Phi}) measures the 
probability, for a randomly chosen closed path of length 2k on our infinite dimensional 
hypercubic lattice, to have area equal to I. The moments of the previous probability 
distribution are defined as 



fe(fe-l)/2 

M« = (/•>*= £ i s P k (l), 

1=0 

and can be computed from the following generating function 

G k (q) 



(23) 



Pk{q) = {q l )y 



Indeed, 



9=1 



G k (l) 







d\ogq 



Pk{q) 



(24) 



(25) 



9=1 



The moments of the deviations with respect to the mean ' (central moments) are 
given by a similar expression: 



m 



dlogq 



Pk(q) q 



-Ml 



(k) 



(26) 



9=1 



The computation of such moments, besides being interesting on its own right, is (at 
least for the first two) a necessary step towards the evaluation of the asymptotic form 
of probability distribution (B^) in the limit k — > oo. We therefore need a systematic 
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procedure to calculate Gk(q), together with its derivatives, at q = 1. Indeed, from 
expression (H), we essentially need to evaluate the derivatives at q = 1 of the integration 



measure (12). Unfortunately such a limit, due to the presence of the Jacobi Theta 
Function, is highly singular: for fixed values of x, v q (x) develops an essential singularity 
at q — 1. A prescription should be given, in our situation the natural one being to chose 
to approach q — 1 from real, smaller than 1, values of q; in a more convenient, equivalent 
notation, we want to evaluate the behaviour of v q (x) in the limit e — > + , e real, defined 
by q = e~ e . With such a choice, as shown in the Appendix, we can write any given 
moment simple combination of trivial Gaussian moments: 



M[ k) 



(2k 



+oo 



dx x 2k c s (x) e 



(27) 



where c s (x) is an even polynomial of degree 4s in x. The only practical difficulty is the 
explicit computation of the s'th coefficient c s (x) of the Taylor expansion around e = + , 
of the function 



(21 



The expansion can be done by hand up to the first few orders in e, but soon becomes 
intractable. It is however a trivial task if we use some symbolic manipulation program 
such as Mathematical^]. We shall give here explicit results for the first few moments: 



M, 



Mj 



(k) 



M, 
M, 



(k) 

2 

(k) 



(29) 



k(k - l)/6 

fc(A;-l)(5A; 2 -A; + 12)/180 
k(k - l)(35/fc 4 + 14A; 3 + 235A; 2 - 188A; + 24)/7560 

k(k - 1)(175A; 6 + 315/fc 5 + 2341 A; 4 - 1959fc 3 + 2056A; 2 - 5664A; - 2160)/226800 

Analogously, the evaluation of the central moments m^ k \ is reduced to the explicit com- 
putation of the s'th coefficient of the Taylor expansion around e = + , of the function 

— « 1 



(k) 



vT 



with the following results: 



-x 2 j e efc(fc-l)/6 



(30) 



in 



o 

(fc 
m\ 

(fc 
(fc 



(fc 



1 


(A; + 3)k(k - l)/45 

(2fc + 3)(2fc + l)A;(fc-l)/945 

k(k - l)(7fc 4 + 37fc 3 + 7k 2 - 108k - 45)/4725 



(31) 
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By direct inspection of the formulae, it is possible to extract the large k behaviour 
of generic moments. Since c s (x) is a polynomial of degree 4s in x, M( fe ) ~ (2k + 4s - 
l)\\/(2k — 1)!! = 0(k 2s ) as k — ► oo|. Following a similar line of thought, and taking 
into account some cancellations which may occur more or less severely, according to the 
parity of s, it is also possible to guess that = 0(k^ s ^ Lp -) for asymptotically large k. 

(k) 

Let us further notice that, on the infinite dimensional hypercubic lattice, from M{ ~ 
k 2 /6, the area of random loops increases on average as the square of their length. This 



should be contrasted to the two-dimensional lattice [ I2| , where the average area increases 
of course as the square of the average linear extension of the loop, that is as the first 
power of the length itself. 

We shall conclude this Section by considering coefficients G k> i for fixed values of I. We 
then have for each integer value I an infinite succession of integer numbers. From direct 
inspection of formula (|i~7|), for asymptotically large values of k we can write: 



k l 

Gk,l — Gkfi -7T 



1 " [ \ 



I fixed. (32) 



6 Asymptotic Behaviour of Distribution Pjfe(Z) 

Let us turn back to probability distribution Pk(l), eq. ( p2]) . We want to investigate its 
asymptotic behaviour as k — > oo. The standard procedure is to rescale the discrete index 
I to a continuous variable t using the average and the standard deviation as computed in 
Section |5| for generic A;: 

^ = M [ k) = M*zI2 ; 
I —> fJ>k + crkt , _ {k) _ {k+ z )k{k -i) (33) 

a k = m 2 — 45 • 



It is now possible to define the asymptotic probability distribution (|22|) as: 



P(t) = lim a k P k (n k + a k t) . (34) 

Unfortunately the coefficients G kj i, are defined in terms of alternate sums of very large 
numbers (about y2l terms, each one of order 0(k k ) } summing up to a much smaller 
number 0), so that the evaluation of the asymptotic distribution P(t) is not completely 
straightforward. Indeed, we are interested into the behaviour of the G k /s for large k and 
large I, I ~ k 2 /6, and the simple estimate of eq. (^) is now completely useless. 

In principle the asymptotic behaviour of coefficients G k (q) for large value of k is 
encoded in the behaviour of generating function R(z,q), eq. (p0|), in the vicinity of its 



"More precisely, it is easy to show that = (fc 2 /6) s [l + 0(l/k)] for large k. 

**Just as an example, let us quote that Gioo,iooo is expressed through formula ( |l7| ) as an alternate 
sum of 45 numbers of order up to 10 200 , being itself of order 10 179 i.e. 21 orders of magnitude smaller. 
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closest (with respect to the origin) singularity in the variable z. The latter is situated 
at z = z c = (1 — q)/4. But since we are in fact interested into the behaviour of G^-,/ 
where I is also asymptotically large, but in a controlled way, see (|33|), we must inspect 
the generating function in the vicinity of its closest (again, with respect to the origin) 
singularity in variable q, that is near q = 1. But in this way z c is collapsing towards the 
origin, conspiring to the construction of an essential singularity in z = z c = as q — ► 1~. 
To tackle this difficulty we shall not investigate the generating function R(z,q), but 
we shall rather resort to the Laplace transform method, applied directly to probability 
distribution Pfc(Z); let us first introduce the Laplace transform 

fc(fc-l)/2 fi-^\ 

Ms) = E P ^) e" V '» ) (35) 

of probability distribution Pk(l), suitably shifted and rescaled according to eq. (|33|). In 
the large k limit, <j>k{s) is expected to converge pointwise to the Laplace transform of 
distribution P(t): 

f+OO 

/ dt P(t) e~ st . (36) 

Since we may write 

_ s 

Ms) =e°* (e 'k) = e v , (37) 

and afc ~ k 3 ^ 2 for large fc, the asymptotic behaviour of Pk(l) is ruled by the behaviour of 
Gk(q) in the neighbourhood of unity. Indeed, after identification of e with s/a^, we can 
make use of the machinery developed in Section || and Appendix A. We start by using 
the fact that, up to terms vanishing exponentially for small e, (q = e~ e ) 

/+oo 
dx x 2h F(x, e) , (38) 
-oo 

with T{x, e) given by equation (|5j). For e = (q = 1), the evaluation of the integral 
reduces to a straightforward application of the saddle-point method: we have to sum up 
the contribution of the neighbourhoods of two stationary points situated in x = x± = 
±y/2k, respectively. We readily get 



G fc (l) = V2 (2k) 



k -k 



(39) 



which is indeed the leading asymptotics for {2k — 1)!!, as it should. Now, for e small but 
positive (we are therefore specializing to positive real values of s), for each of the two 
stationary point x — x±, we just shift the integration variable x to x±+y and expand the 
integrand, including !F(x,e), around k = oo. Considering e.g. the contribution coming 
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from the positive stationary point x + , and taking into account only terms up to 0(\ogk), 
we get 

V27T J 

Of course the situation in which s < has to be computed, too, but the result turns out 
to be independent of the sign of s. Integrating, and adding the analogous contribution 
coming from the region in the neighbourhood of x_ , we finally obtain 



G k (e't) = v/2 (2k) k e~ k e^"^ 
and thus, for asymptotically large values of k 



l + Ol 



(41) 



(f>k{s) 



e 2 



l + 0( 



(42) 



Indeed, by comparison with eq. (|3"E|), we have thus proven that P(t) as defined in equation 
341) is a standard Gaussian distribution: 



Pit) 



'Ik 



(43) 



This result gives us further information about the asymptotic behaviour of the central 
moments of distribution P&(Z): indeed, rescaling each central moment my^ by a 
factor a^ s ~ k~z s (correspondingly to the rescaling of the discrete variable I to the 
continuous one t), and performing the k — > oo limit, we must find the moments of a 
standard (unit variance) Gaussian distribution, namely or (s — 1)!!, for odd or even 
values of s, respectively. This on one hand confirms rigorously the estimate = 
(T^ s ■ (s — 1)!! [1 + 0(l/k)] (s even) and on the other hand matches the guessed behaviour 
m (*0 k^ s ^- p - (s odd). 

As for the area versus length distribution, the result we have found, eq. 



does 



not agree with the behaviour proposed in Refs. |IT], [12] ; however the role of the D — » oo 
limit, inherent to our treatment, must be taken into proper account. In this limit, all 
loops with more than one plaquette of the subtended surface lying over the same plane 
are neglected; on the other hand, on the two-dimensional lattice considered in Refs. 
11, |12| there is only one possible plane. Moreover, having averaged over orientations, we 



restrict ourselves to the absolute value of area, while in |TT|, |T2j , algebraic (oriented) area 
is considered. 

When considering very compact loops, their distribution is ruled by the left tail of 
eq. (|4"3T). In Ref. |H| the different but somehow related problem of the distribution 
p(s, t) of random walks of t steps visiting exactly s different sites was addressed on a 
cubic lattice. In particular, an asymptotic formula in the regime of large t and relatively 
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small s was proposed, which describes the distribution of compact walks. However the 
two distributions are not comparable: indeed, in our situation of infinite dimensionality, 
the most compact loops of length 2k have zero enclosed area, but still visit a relatively 
large number of distinct sites, at least k + 1; more precisely, when the dimensionality 
of the lattice is very large, the number of walks which take more than one step in the 
same direction is of order 1/D with respect to those stepping in k different directions, 
and therefore even the most compact walk, the one coming back to the origin every two 
steps, at leading order in D visits k different sites, in addition to the origin. In other 
words, the connection between the two notions of compactness, as small enclosed area, 
or as small number of visited sites, is lost in the limit of infinite dimensionality. 



7 An Infinite Hierarchy of Nontrivial Identities 

We would like to conclude with a simple consideration, which in our opinion is however 
rather striking and deserve further investigation. The exact evaluation of moments 
gives as a byproduct an infinite set of nontrivial and somehow intriguing identities, the 
simplest of them beeing: 

fc(fc-l)/2 

J2 G kil = (2k-1)\\ k = 0,1,2,... (44) 



and in general 



A:(fc-l)/2 

£ G Kl T = M^ (2k-l)\\ k = 0,1,2,... (45) 



with Gk,i and explicitely given by formulae (|17]) and ( p7| - p8|) , respectively. These 
identities are of course encoded in the Jacobi Theta function, essentially in its modular 
invariance properties, eq. (|4*8|), but they come nevertheless rather unexpected. 

Other sets of identities can be obtained if we are somehow able to compute the Gk,i in 
a indipendent and simple way. This is indeed very easy e.g. for Gk,k(k-i)/2> which counts 
the number of diagrams with 2k external points and maximal crossing number, and is 
thus obviously equal to 1 for any k (just connect the j th point with the (j + k) th one, 
j — 1, . . . , k). From eq. ([IT]) we may therefore write: 



=0 



(46) 

Identities of this second sort are however more obvious: eq. (|4~6|) is just the simplest 
example of a whole bunch of identities encoded in an obvious simmetry of polynomial 
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Gk(q) : by applying the Theorem of Section |3| (eq. QI6D ) to the new polynomial Gp.(q) = 
q k ( k ~ x '' 2 Gk{l/q) } with Gk(q) expressed in terms of eq. ([13]) one obtains new coefficients 
Gk,i, which obviously satisfy G/~,i = Gk,(k(k-i)/2)-u setting in particular I = 0, we readily 
get eq. ©. 



8 Conclusions 

The problem of counting closed paths of given length on a lattice, according to the 
area they enclose, is a difficult one. In the present paper we have tackled the question 
on an hypercubic lattice, in the limit of infinite dimensionality. Moreover an average 
over orientations is inherent to the method we have used. But these conditions have 
allowed us to give a closed answer, which is moreover exact under the restriction we 
have assumed. The main flaw of the whole method resides in our opinion into the fact 
that the computation of even only the first corrections in 1/D to coefficients Gk,i is out 
of reach. On the other hand, we have presented an explicit solution to a combinatorial 
problem which could hardly be addressed directly; in this respect, the connection between 



the abstract combinatorial puzzle and the physical model proposed in [13j is crucial. 
The exact solution proposed for coefficients Gk,i has been investigated in detail; their 
asymptotic behaviour has been evaluated; finally, as a byproduct, an infinite set of new 
nontrivial identities has been obtained. 
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Appendix A 

We want to build up a procedure to compute in a systematic way derivatives of 
arbitrary order of Gk{q) with respect to logg, at q — 1. From the explicit expression ( |14D 
of Gk{q), we have to take into account both the q dependence of the integration interval 
and of the integration measure v q {x), given by eq. (|T2|). As for the latter, the task 
is indeed non completely straightforward, since for fixed values of x, v q (x) develops an 
essential singularity in the limit q — ► 1, and should somehow be regularized. In presence 
of an essential singularity, the result of the limiting procedure depends from the direction 
chosen to approach the singularity, and a regularization is a prescription on how the 
singularity is approached. Indeed, if we restrict q to the real axis and approach q = 1 
from smaller values, it is possible to write the integration measure u q (x) (x = -e= cos 6) 
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as a sum of a completely regular contribution, Taylor expandable in the real variable 
e, (q = e~ e ) plus a singular contribution vanishing with all its derivatives for e — > + . 
Indeed we shall show that: 



Vq (x) = J 7 {x,e) \l + 0(e~ 1/e ) 



(47) 



T{x, e) being an analytic function of e. 
Using modular invariance, we can write 



1 



\7T / \J—IT \7TT T 

on the other hand, from definition (|T3D, letting r = j-, 



(48) 



6i 



-1 




1 


V 7TT 




i 



e-ire^ £ (-1)" e^"'^ 5 **) 



(49) 



neZ 



For arbitrary small positive e, and for - in the open interval (0, 1), the leading contribution 
is given by the n = term in the infinite sum, all other terms being exponentially small. 



ei 



for r 


te 

~ 2tt 


, and < e < 1 


( £ ' 




1 7T 2 2-n-e 

~ — 6 2e e e 


\ TIT 




' % 



l + Ofe « 



(50) 



and for the measure ^ g (x) we readily get, in the limit q — *■ 1 , or g = e £ , 6 ^ + , the 
approximate behaviour (f47|), with 



e~ £ es 



v/2 



7T6 



vT 



e~ e es 



v/2 



7T6 



exp 



exp 



e V 2 



arcsm" I — yl — e 
e V2 



^7T 



e 2 



E 

r=0 



Cr(a?) P 



'/"' ! 



(51) 
(52) 
(53) 



where the coefficients c r (x) are given by the Taylor expansion in the neighbourhood of 
the origin e = of the analytical function: 

1 



Vl — e~ e es — e 2 > 



»( f> /l^-x»] 



(54) 



Co(x) = 1 trivially while higher coefficients c r (x) can be shown to be even polynomials 
of degree 4r in x. In the neighbourhood of q = 1 we may therefore write, up to terms 
vanishing exponentially for small e: 



G k (q) 



'2tT J- 



/ ' 2 2 00 

^ dx x 2k ■ 



c r (^c) , r 



(55) 
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Let us now observe that the integration limits diverge as 1/ y/e, while the integrand falls 
down as a quadratic exponential. Therefore, up to terms which vanish exponentially for 
small e, the integration interval can readily be extended to the whole real axis, and the 



computation of moments ( Pq) is finally reduced to the Taylor expansion of ([54]) and the 
evaluation of a sum of standard Gaussian integrals: 

/ "1 \S j «_J_QQ 2 



As an example, let us compute explicitly the first moment, M\ . The first coefficient 
of the Taylor expansion of (B5) is 



Ci(x) = 1 — x 2 x 4 , (57) 

u > 8 4 24 ' v ; 

and therefore, denoting with ((. . .)) the Gaussian average (with variance equal to 1), and 
exploiting the standard relation ((x 2k )) = (2k — 1)!! 



Mx(k) 



8 4 V 
k{k- 1 



l((x 2k ))-\((x 2k+2 )) + ^((x 2 ^)) 



(5? 



l-\(2k + l) + ^(2k + 3)(2k + l) (59) 

(60) 
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